Indirect Dissociative Recombination of LiHt + e 
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^Department of Physics and JILA, University of Colorado, Boulder Colorado 80309 

We present the results of calculations determining the cross sections for indirect dissociative 
recombination of LiHj + e~. These calculations employ multichannel quantum defect theory and 
Fano's rovibrational frame transformation technique to obtain the indirect DR cross section in 
the manner described by Ref.[3|. We use ah initio electron- molecule scattering codes to calculate 
quantum defects. In contrast to Hj, the LiHj molecule exhibits considerable mixing between 
rotation and vibration; however, by incorporating an exact treatment of the rovibrational dynamics 
of the LiHj", we show that this mixing has only a small effect on the observed DR rate. We calculate 



a large DR rate for this cation, 4.0 x 10 cm s 
PACS numbers: 03.65.Nk, 34.80.-i, 34.80.Lx, 33.20.Wr 

I. INTRODUCTION 

Dissociative recombination, the process by which a 
cation recombines with a free electron and dissociates, 



AB^ 



A + B 



(1) 



has received much theoretical and experimental interest 
in the past two decades [3, y] . Innovations at both sides 
of the scientific process have spurred this interest. The 
development of storage ring experiments [J] has been the 
key innovation on the experimental side. Storage rings 
allow the preparation of cation species that are rovibra- 
tionally cold, such that a small number of initial rovi- 
brational states are populated. Such devices also enable 
the synchronization of cation and electron beams, such 
that the relative kinetic energy between the two can be 
precisely controlled. As a result, DR rate coefficients 
can be determined with unprecedented resolution, and 
structures in the rate coefficient as a function of relative 
kinetic energy may be elucidated. 

The current theoretical understanding of the dissocia- 
tive recombination process provides two mechanisms by 
which it may occur. These mechanisms are labeled the 
"direct" and the "indirect" process. The direct pro- 
cess involves temporary capture of the electron into a 
metastable electronic state of the neutral. Such resonant 
electron capture was pointed out by Bates in lOSOQ^and 
quantitatively formulated later by O'Malley in 1966[g]; it 
is particularly effective in capturing low-energy (thermal) 
electrons when the Born-Oppenheimer potential energy 
curve of the metastable neutral state crosses the curve 
of the ground state cation species within the Franck- 
Condon region of the latter. It may also be the only 
viable mechanism of dissociative recombination at high 
incident electron energy. 

When there is no Born-Oppenheimer curve of the neu- 
tral that crosses within the Franck-Condon region of the 
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at 1 meV incident electron energy. 



cation, it is the indirect mechanismQ that is responsible 
for any observed dissociative recombination. The indi- 
rect mechanism is favored by low kinetic energy of the 
electron-cation collision. The indirect process, like the 
direct process, is a resonant phenomenon; however, in 
this case the resonances arc rovibrational Feshbach reso- 
nances, not electronic resonances as in the direct process. 

Until recently, the consistent, accurate mathematical 
and numerical description of the indirect mechanism was 
elusive Q. Perhaps the most vexing problem was that of 
the dissociative recombination of H;^ , because the disso- 
ciative recombination of this species plays an important 
role in interstellar chemistry, and due to the numerous 
failures of theory to accurately predict the rate observed 
by experiment. Adding to the mystery was the consider- 
able spread in experimental results, rangingfrom 2.3 x 
lO"'^ to less than IQ-^" cm'^ g-i at 300°Kg. 

However, the theory outlined in Ref.[l|, involving a 
frame transformation with Sicgert states representing the 
outgoing dissociative flux, has been applied to several 
systems and has thus far shown consistently good results 
in predictin g in direct DR rates. A series of theoretical 

works [3, [l3i llHj ll2l O'^ t^^ ^^ of ^t ^^^ isotopomers 
obtained unprecedented agreement with experiment for 
this difficult system, matching both the overall magni- 
tude and most of the structure of the experimental cross 
section[lj, [ij, [l^. Further use of the method has in- 
cluded a study of LiH+ [1^, \VR i: hat reproduced the ex- 
periment of S. Krohn et al. [T^ \lM extremely well in 
all but the lowest part of the measured incident energy 
range. 

In the present article, we examine the dissociative re- 
combination of another species, namely LiHj . Despite 
any superficial similarity to H;^, the two cations are in 
fact quite different, and we view these calculations as a 
further step toward validating and generalizing the the- 
ory. In particular, the rovibrational structure of the LiHj 
cation is more complicated than that of H;^; whereas 
Rcf. 9], and later, Ref.[l2| obtained excellent agreement 
with experiment by using a rigid rotor approximation 
for the vibrational states of H;^, the LiH^ cation is 
well described as a Li"*" cation weakly bound to an II2 
molecule, which fragments may rotate relatively indepcn- 



dently. Thus, in the present work we incorporate an ex- 
act treatment of the rovibrational Hamiltonian and com- 
pare it to a rigid- rotor treatment. This work represents 
the first such exact treatment of the ionic rovibrational 
motion for indirect dissociative recombination in a poly- 
atomic species. 

This paper is organized as follows. We briefly intro- 
duce the electronic structure of LiH^ and LiH2 in Sec- 
tion mi W e use the the Swedish-Molecule and UK R- 
matrix[20| codes to calculate fixed-nuclei electron scat- 
tering quantum defect matrices, and we describe these 
calculations and present the results in Section IIIIl A 
description of the calculation of the rovibrational states 
of the cation, including an explanation of the coordinate 
system we use, comprises Section|TVl In Section|V]we de- 
scribe how we account for the outgoing dissociative flux; 
we employ a method different from that used in previous 
calculations, using exterior complex scaling 2l| instead of 
Siegert states to enforce outgoing wave boundary condi- 
tions on the vibrational basis. In Section IVll we describe 
the rovibrational frame transformation and explain how 
the nuclear statistics are taken into account. Finally, in 
Section FVlII we present the calculated cross sections. 



II. ELECTRONIC STRUCTURE OF LIHt AND 



LIH2 



The ground electronic state of the LiH2 molecule and 
the LiH^ cation are well described qualitatively as a Li 
atom or Li"*" cation weakly bound to an H2 molecule. 
Both states have an equilibrium geometry with equal Li- 
H bond lengths, and in such a geometry the molecule 
belongs to the C2t, point group. Using the labels appro- 
priate to C2i, symmetry, the electronic configuration of 
the cation is la^ 2a^, for overall ^Ai symmetry. The addi- 
tional electron for LiH2 goes into the 3ai orbital (approx- 
imately the Li 2s orbital). When the Li-H bond lengths 
are unequal, the molecule belongs to the Cs point group 
and the cation configuration is labeled la'^ 2a'^. 



Prior calculations [23, |23|, M, M, M, M, M, ]2M have es- 
tablished that the equilibrium geometry of the cation has 
rHH = 1.42ao and R = 3.62ao, where R is the distance 
between the Li and the H2 center of mass. The two body 
asymptote Li+ + H2 lies only 0.286eV higher [s^. The 
three body asymptote (Li+ + H + H) hes much higher, 
5.034eV [sol- The LiH-l- complex is weakly bound with 
a dissociation energy of 0.112eV and therefore this two- 
body breakup channel is essentially isoenergetic with the 
three-body channel. 

Because the excitation energy of the Li+ cation is very 
high ~ 60.92eV - the lowest-lying electronic excitations 
of LiH^ correspond to states of the Li atom bound to a 
H^ molecule. The ionization energies of Li and II2 are 
5.39 and 15.43eV, respectively, and so we expect the first 
excited state to occur at roughly lOeV. 
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FIG. 1: Quantum defect matrix elements (top two panels) 
and cation potential energy surface [301] (bottom panel) as a 
function of the Jacobi coordinate R, fixing run = 1.4ao and 
7 = 90°. Diagonal matrix elements are labeled with a single 
partial wave; coupling matrix elements are labeled with the 
two. 



III. FIXED NUCLEI SCATTERING 
CALCULATIONS ON e + LIH+ 

The first step in the present treatment is the calcula- 
tion of the fixed-nuclei quantum defect matrices, in the 
body frame, which describe the scattering of an electron 
from the LiHj cation, with the positions of the nuclei 
frozen in space. To perform this task we employ the 
polyatomic UK R-matrix scattering codes [201 based on 
the Swedish-Molecule electronic structure suite. 

The R-matrix calculation is defined as follows. We 



employed an augVTZ STO basis set[3lj and a 20 bohr 
spherical R-matrix box radius. The center of mass of the 
Lill2 cation was placed at the origin. We first perform 
a Hartree-Fock calculation on the cation using the la'^ 
2a'^ configuration. The target wavefunctions are defined 
as having the la' orbital (the Li Is orbital) frozen in 
double occupation, with the remaining two electrons dis- 
tributed among the space 2-6a' and la". We keep the first 
nine roots of this complete active space configuration- 
interaction (CAS-CI) calculation to include in the scat- 
tering calculation. These correspond to the ground state. 
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FIG. 2: Quantum defect matrix elements (top two panels) 
and cation potential energy surface [301] (bottom panel) as a 
function of the Jacobi coordinate thh, fixing R = 3.62ao and 
7 = 90°. 



FIG. 3: Quantum defect matrix elements (top two panels) 
and cation potential energy surface [30] (bottom panel) as a 
function of the Jacobi coordinate 7, fixing R = 3.62ao and 
THH — 1.4ao. 



and excited states that correspond roughly to an Hj 
molecule bound to a Li atom in its X ^5 or ^P con- 
figurations, singlet or triplet coupled. Thus we have four 
^A' states, three ^A' states, and one ^A" and '^A" state. 
At the equilibrium geometry of the cation our treatment 
places these states between 12.86 and 17.13eV. 

To the target orbital space we add a set of uncontracted 
Gaussians that represent the scattering electron. This set 
is obtained using the UK R-matrix code GTOBAS[33]. 
which optimizes the set to best fit a set of coulomb wave- 
functions orthonormal over the R-matrix sphere. We in- 
clude 15 s orbitals. 13 p orbitals, and 12 d orbitals opti- 
mized to fit coulomb wavefunctions up to 10 hartree. 

The five-electron space included in the R-matrix calcu- 
lation is defined as follows. We include a close-coupling 
expansion corresponding to the first nine states discussed 
above times scattering orbitals, plus penetration terms in 
which all five electrons are distributed among the target 
orbitals. again keeping the la' orbital doubly occupied. 
The calculation is performed in overall A' or A" symme- 
try. 

These calculations yield the fixed-nuclei quantum de- 
fect matrices fJ.im,i'm' that are included in the later steps 
of the dissociative recombination calculation. The quan- 



tum defect matrix is defined in terms of the fixed-nuclei 
S-matrix as /i = — 2^ ln(5). These quantum defect ma- 
trices depend weakly on the incident electron energy; 
we evaluate them at an incident electron energy of two 
nieV. We construct an interpolated quantum defect ma- 
trix by splining the calculated quantum defect matrices 
over the Jacobi coordinate range 2.4ao < R < 5.6ao, 
0.6ao < thh < 2.4ao, all 7. 

Plots of the splined quantum defect surfaces are shown 
in Figs. [U [21 and[3] These figures show three cuts through 
the quantum defect surfaces and the corresponding cuts 
through the cation potential energy surface; all three 
points contain the point (_R=3.62ao, rnH—^-^o-o-, 7=90°) 
in Jacobi coordinates. The cuts are in the R direction 
(Fig. [1]), the Thh direction (Fig. [2]), and the 7 direction 
(Fig. [2|). The convention in these figures is that all of the 
diagonal quantum defects are labeled and labeled with a 
single channel index, and some of the off-diagonal defects 
are labeled and labeled by the corresponding pair of in- 
dices. The molecule lies in the yz plane and the vector 
R, which connects the Li atom to the H2 center of mass, 
is coUinear with the z axis. 

For the calculation in overall A" symmetry there are 
three electronic channels included in the R-matrix calcu- 




FIG. 4: "R-embedding" rovibrational Jacobi coordinate sys- 
tem with origin at the center of mass. Primed and unprimed 
axes refer to BF and SF frames, respectively. The BF x' z' 
and x'y planes are both marked with a thin line circle and 
the SF xz and xy planes are marked with dashed circles. The 
line of nodes is also drawn. The molecule resides in the BF 
x' z' plane. The Delves hyperspherical coordinates Q and ,^ 
used in the DR calculation are defined in terms of R and r. 



lation, the p^^ d^^, and dxy We find that the quantum 
defects in A" symmetry are relatively small. For the cal- 
culation in overall A' symmetry there are six electronic 
channels included in the R-matrix calculation. The quan- 
tum defect matrix elements involving pj, and d^^-y^ are 
small relative to the other four. 



IV. CALCULATION OF BOUND AND 

OUTGOING WAVE ROVIBRATIONAL STATES 

OF THE CATION 

The next step in the DR treatment involves the cal- 
culation of rovibrational eigenfunctions using the ground 
cation potential energy surface. We employ the surface 
of Martinazzo et al. [30| , which includes the proper long- 
range behavior of the potential. 



A. Coordinate system and Hamiltonian 

As in previous treatments[3, [lO, E, [23 1 ^^ ^^e a hy- 
perspherical coordinate system and construct rovibra- 
tional states in an adiabatic hyperspherical basis |34l|. 
The adiabatic expansion helps to reduce the size of the 
calculation. 

In contrast to the previous treatments we use Delves 
hyperspherical coordinates [3a, IM]. These coordinates 



are built from the Jacobi coordinate system appropri- 
ate to the system, in which r^H denotes the H2 bond 
length, R denotes the distance between the Li atom and 
the H2 center of mass, and 7 denotes the angle between 
the two corresponding vectors. The Delves coordinates 
consist of the Jacobi coordinate 7, plus two additional 
coordinates ,'% and 6, 



9 — tan 




Mr r 

For calculations with nonzero total cation rotational 
angular momentum J+ , we employ the i?-embedding co- 
ordinate systeml37| in which the Euler angles a, f3, C ori- 
ent the molecular z' axis, coUinear with the R vector, and 
the molecular x'z' plane, which contains the molecule, 
relative to space-fixed axes. This coordinate system is 
depicted in Figure [H 

We employ the exact rovibrational Hamiltonian for this 
coordinate system, taken from its form in Jacobi coordi- 
nate system - see, for example, Refs.j33l.l38j. 
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(3) 



In this equation, the operators p and j± arc the to- 
tal and raising/lowering operators of the diatom angular 
momentum. This Hamiltonian operates on the expansion 
coefficients xk in the following expansion of a wavefunc- 
tion. 



*j+A/ = Y^ 



K 






D 



MK 



ia,PX) , 



(4) 



where the basis of Dfjj^la, (3, () is the set of normalized 
Wigner rotation matrices (and BF angular momentum 
eigenstates) 



^M/^(",AC) = / 



^"^^^mV(«,/?,c) 



(5) 



B. Coupled adiabatic hyperspherical treatment 

The first step in calculating the rovibrational states is 
to calculate the adiabatic hyperspherical basis. There- 
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FIG. 5: Rovibrational energy eigenvalues for J^=0 and 1 
(relative to three-body breakup) calculated presently with 
the surface of Martinazzo et al. [Sfl and the full rovibrational 
Hamiltonian, Eq.Q (filled dots); those calculated for J'^ — 1 
with a rigid rotor approximation (empty dots); and those cal- 
culated by Sanz et a/. [4^ with the Martinazzo surface (trian- 
gles). 



where the matrix H is defined 



ttJ+ ^ ,J+ 



i^^) + 04.'/ (7k). 



(10) 



where T^ is the Gauss-Lobatto kinetic energy matrix for 
tlie hyperradius, and where the matrix O"^ is the overlap 
matrix 



or 



W' = (xf (^■0|x/'^(^.')) , (11) 



brackets denoting integration over aU degrees of freedom 
except M. 



Rigid rotor approximation and rovibrational 
energies 



fore, defining H'^^ = Hf{^) - ^-5/2_i_^^5/2 
where Hq is the adiabatic Hamiltonian, we first solve 
for adiabatic basis functions Xj *''(^i7i ctj/^j C;-^) a-nd 
eigenvalues e^ (M) , 



'mxf{0,j,a,i3,C;m, 



(6) 



where we expand x 



J+M 



as 



xf ^(0,7,«,/3,C;^) = $3 



K 



^5/2 



(7) 

The a-th rovibrational eigenfunction for total cation 
rotational angular momentum J"*" is then expanded as 

^j+Ma = E 4^o.Mmxf"i9, 7, «, /?, C; ^0 ■ (8) 



The coefficients c/^-, multiply basis functions (f>i{^) based 
on gridpoints Mi. These functions comprise a Discrete 
Variable Representation (DVR) [33, So, Ej , specifically, 
the Gauss-Lobatto finite element DVR[44| with five el- 
ements 1.6 bohr long, starting at 2.0 bohr, and order 
10 within each element. For the hyperangular degree of 
freedom 6 we also use Gauss-Lobatto DVR, but with one 
element, and 60th order. The wavefunction is defined to 
be zero at 6* = and 90°. For the 7 degree of freedom 
we use Legendre DVR based upon associated Legendre 
functions Pik- The potential is evaluated using the DVR 
approximation, which corresponds to a diagonal repre- 
sentation. 

To calculate the full vibrational wavefunctions includ- 
ing the nonadiabatic coupling, we employ the slow vari- 
able discretization of Tolstikhin[44|. and therefore solve 
the matrix equation for the coefficients c,"^ 



2ja ' 



H'^'c^ 



Ei^cr 



(9) 



To calculate the rigid rotor states, we calculate the 
vibrational states for total cation angular momentum 
J^ = 0, obtaining their wavefunctions ^ooa ^^^d ener- 
gies E^. We find the principal moments of inertia A, B, 
and C for each state; the largest of these, A, is perpendic- 
ular to the molecular plane. We use this moment as the 
axis of quantization and then diagonalize the asymmetric 
top hamiltonian 



H, 



B + C 



rigid 



r 



2A-B-C 



4 
B-C 

8 



-Ji 



J+.P 



J-J- 



(12) 



in the basis D'^^j^ for a given total cation angular momen- 
tum J^ . (In this equation, J are raising and lowering 
operators of the projection, K, of the total angular mo- 
mentum on the body-fixed axis of quantization. They are 
not to be confused with the total cation angular momen- 
tum J+, where J+( J+ + 1) is the eigenvalue of the total 
angular momentum squared operator J^. K is the eigen- 
value of Jz- M is, yet again, arbitrary.) For each value 
of J+ and each J+ = state '^QQa^ ^^ obtain 2J+ + 1 
eigenvalues which are added to E^ to yield the rigid ro- 
tor energies for that vibrational state. For the purposes 
of the rotational frame transformation, we transform the 
eigenvectors of Hrigid such that their axis of quantiza- 
tion, conjugate to the eigenvalue K, is parallel with the 
Jacobi vector R, not perpendicular to the plane. 

The J"*" = vibrational energies (which are eqiuvalent 
in the rigid rotor and full rovibrational calculations) are 
in good agreement with the results of Sanz et al. [43| . For 
J^ > 0; The rigid rotor approximation gives significantly 
different low-lying eigenvalues than the full rovibrational 
calculation. In Figure [5] we plot the energies for rovi- 
brational states with J+=0 and 1. The eigenvalues of 
Sanz et al. for J+=0 agree reasonably well with ours. 
For J'^ = l we plot eigenvalues calculated with the full 
Hamiltonian, Eq.(I3]), as well as those calculated in the 
rigid rotor approximation. One can clearly see that it is 
not accurate to treat this molecule as a rigid rotor. 
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FIG. 6: Boltzniann weights at 300° K binned by total cation 
angular momentum J+. 



We plot the Boltzniann weights binned by cation rovi- 
brational angular momentum value J+ in Figure [6l The 
number of rovibrational states goes as (2J+ + 1)^ and 
thus the most probable J~^ value at 300° K is six. 



Siegert states. We will present a more thorough compari- 
son of these different methods of enforcing outgoing- wave 
boundary conditions in a frame transformation calcula- 
tion in a forthcoming publication. 

To calculate the ECS eigenvectors, the final finite ele- 
ment in the ^ degree of freedom is scaled according to 
.^ ^ ^0 + e*^(^ - ^o), where ,'%o is the boundary be- 
tween the fourth and fifth elements at .'% — 8.4ao. We 
employ a scaling angle of ^tt. As with Siegert states, this 
leads to a discretized representation of the dissociative 
Li+ + H2 vibrational continuum in which the outgoing 
wave states have a negative imaginary component to their 
energy. 

Because the coordinate ^ is scaled into the complex 
plane, it is ideal (but often not necessary [43|) to analyt- 
ically continue the potential energy surface V{^,9,'^). 
We do so by ensuring that the long-range components to 
the Martinazzo et al. surface are evaluated for complex 
arguments. We evaluate their switching formula (third 
equation on page 11245 of their publication [SOj') by tak- 
ing the absolute value of the argument. 



D. Nuclear statistics 

The full rovibrational Hamiltonian is invariant with re- 
spect to permutations of the two hydrogen atoms. There- 
fore, the rovibrational eigenfunctions will have an eigen- 
value of either -f 1 or -1 with respect to this permuta- 
tion operation, which can be expressed ( 7 —> 90° — 7; 
C — > C + 180° ). Given that the hydrogen atom is a 
fermion, the -1-1 states are paired with a singlet (para) 
nuclear spin wavefunction, and the -1 states are paired 
with a triplet (ortho) nuclear spin wavefunction. This 
gives the -t-1 and -1 states statistical weights of 1 and 3, 
respectively. 

The full rotational/rovibrational frame transforma- 
tion, described later, does not affect the nuclear statis- 
tics. However, the rovibration-only frame transformation 
mixes states with different permutation eigenvalues, and 
therefore we cannot account for the proper nuclear statis- 
tics with this transformation. 



V. REPRESENTATION OF OUTGOING FLUX 

The previous implementations of the present theory 
have employed Siegert pseudostate s 1451 or complex ab- 
sorbing potentials (CAPs) [4y, [43, [3 to represent the 
outgoing flux corresponding to dissociative recombina- 
tion. In contrast, in the current implementation we 
employ exterior complex scahng (ECS)[2l|, |4^, [50, [5l|, 
[52l . [53I . [54| to enforce outgoing-wave boundary condi- 
tions. We have found that the use of ECS or CAP 
states within a MQDT frame-transformation calculation 
is more straightforward than the use of Siegert states, as 
the completeness and orthogonality relationships of the 
former types of eigenvectors are simpler than those of 



VI. ROVIBRATIONAL FRAME 
TRANSFORMATION 

A. Introduction 

The rovibrational frame transformation comprises the 
central part of the present calculation. Frame transfor- 
mation techniques were originally developed by Fano[55l 
[56[ and have found much use in atomic and molecular 
theory. The central idea of a frame transformation is to 
take an S-matrix, which is labeled by incoming and out- 
going channel indices, and transform that S-matrix to a 
new channel basis. In its simplest incarnation, adopted 
here, this transformation is exact if the fixed-nuclei quan- 
tum defects are constant with respect to energy. The 
transformation is accomplished via a unitary matrix that 
relates the first set of channels to the second. Usually, 
the first set of channel indices are appropriate to describe 
the system when the scattered electron is near the atomic 
or molecular target, and the second set of channel indices 
are appropriate when the electron has escaped far from 
the target. The coefhents of the original rotational frame 
transformation for a diatomic molecule [55| are simply 
Clebsch-Gordan coefficients. Other unitary transforma- 
tions may be applied for different physical situations: for 
the calculation of Stark states [531, to transform between 
LS and J J couplinglSa. [59l. or to transform between 
molecular Hund's cases [GOf. 

The frame transformation is applied to molecular vi- 
bration in much the same way it is applied to rotation. 
When the scattered electron is close to the molecule, it is 
moving very fast compared to the molecular framework, 
and therefore the scattering may be calculated by fix- 
ing the nuclei and obtaining fixed-nuclei, body-frame S- 
matrices sim^i'm'io) where (fare the internal coordinates 



of the molecule and Im, I'm' label the partial wave elec- 
tron scattering channels in the body frame. The frame 
transformation provides that the full S-matrix, which has 
vibrational channel indices as well as electronic channel 
indices, is found via Sima,i'm'i3 = {Xa\sim,i'm'\xi3) where 
the brackets denote integration over the internal degrees 
of freedom q. 

It is important to note that this vibrational 
frame transformation is different from the Chase 
approximation 6l| . The frame transformation is applied 
to the "short-range" S-matrices of Multichannel Quan- 
tum Defect Theory (MQDT) [H, [H, [H , which have in- 
dices including not only open but also closed channels. 
As a result, complicated nonadiabatic effects caused by 
the long-range potential (here a coulomb potential) may 
be accounted for by the theorv|65ll66j. 

The most accurate versions of the vibrational frame 
transformation theory[63, [6a, [£3 incorporate the energy 
dependence of the fixed-nuclei S-matrix. We do not do so 
and instead evaluate the fixed-nuclei S-matrices at 2meV, 
implicitly making the assumption that these S-matrices 
are constant with respect to incident electron energy. 

We note that many other treatments of dissociative 
recombination within MQDT have been devised. See, for 
example, Refs.[70l.l7H. However, the current formulation 
is perhaps the most easily applicable to a polyatomic 
molecule. 



B. Rovibrational frame transformations for the 
asymmetric top 

Child and Jungen[72| have already derived the rota- 
tional frame transformation for the asymmetric top. We 
perform both a rovibration-only frame transformation 
and a rovibrational/rotational frame transformation that 
uses their result. 

For the vibration-only frame transformation we calcu- 
late 






W- 



J+Ma\ ^lm,l'm' 



^J+MI3 



(13) 



where value of the index M is irrelevant. 

The rovibrational frame transformation of Child and 
Jungen[73| will not be repeated in full detail here. It com- 
prises a square unitary transformation matrix for each 
value of J (total angular momentum) and / (the angu- 
lar momentum of the electron). It transforms from the 
body-fixed representation, with quantum numbers m and 
K - denoting the projection of the electron angular mo- 
mentum about the molecular axis and the projection of 
total angular momentum - to the space-fixed representa- 
tion, with quantum numbers J"*" and X+, denoting the 
total angular momentum of the cation and its projection. 
The body-fixed S-matrices are independent of K. Thus, 



^U+K+,l'.J+'K+'\^^^^n — Z^ UmK,J+K+ 

mKm' 



The full rovibrationally and rotationally transformed S-matrix is then 
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The index M is again irrelevant. 
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C. Channel closing and dissociative recombination 
cross section 

The final step in the present theory is the construc- 
tion of the physical, open-channel S-matrix in terms of 
the closed-channel S-matrices calculated from the frame 
transformation. Whereas the latter are assumed to be 
energy-independent, a strong energy dependence is in- 
troduced to the former by the formulajOf 



<--^(-C'j — ^oo ^oc \Jcc 



'2il3{E) 



Sr_ 



I3{E) 



nSi. 



(16) 



we introduce the notation S^ for the physical S-matrix. 

Because the higher-energy rovibrational states lie 
above the dissociation energy to Li+ + H2, they have 
outgoing-wave components and negative imaginary com- 
ponents to their energy. As a result, the physical S- 
matrix is subunitary and we assign the missing part to 
dissociative recombination. Thus, for the vibration-only 
transform, we sum over the contributions of each partial 
wave in the electronic channel. 



I'm'P 



Ira 



, (17) 



where the subscript c and o denote the closed and open 

channel subblocks of the MQDT S-matrix S'^ or S'^ , and and for the full rotational plus vibrational frame trans- 



formation, 






J+laJ+'l'P 



(18) 



where a and J+ denote the initial rovibrational state. 

We Bohzmann-avcrage these results, assuming a 
cation temperature of 300° K. Thus Q, 



a,^biE) = ^ E (2-/+ + 1) ^a^ iE)e^^ 



J+a 



O'rot 



(^)-^E 



JJ+a 



2J+1 
2J^ 



Y '^j^jEy 



tl±o 



E=J2 (2J+ + 1) e^i^ 



(19) 



(20) 



(21) 



J+c 



with T=300° K. 

Finally, we convolute the results with respect to the un- 
certainty in the incident electron kinetic energy. For the 
present results we use a standard deviation of \/2 meV in 
both the parallel and transverse directions, and perform 
the averaging as described in Ref . |17| . 



VII. RESULTS: DISSOCIATIVE 
RECOMBINATION CROSS SECTIONS 

We seek to determine how relevant the inclusion of 
the exact cation rovibrational dynamics is to the exper- 
imentally observed DR rate. The raw DR cross sections 
that we calculate show considerable structure that de- 
pends upon whether an exact or rigid-rotor treatment 
of the rovibrational dynamics is used. However, experi- 
ments operate with a thermal sample of cation targets, 
including many rovibrational states, and use a beam 
of electrons with a small spread in energies. Storage- 
ring experiments are performed with cool cation targets, 
with rovibrational temperatures typically on the order 
of 300° K. In order to compare with results obtained 
under these conditions, we Boltzmann-average over ap- 
proximately 300 initial rovibrational states of the LiHj 
cation, and account for the uncertainty in the incident 
electron energy, taken here to be 2meV (v^nieV in the 
parallel and transverse directions). In doing so, much of 
the structure in the DR cross section is lost, and we find 
that the rigid rotor treatment is probably sufficient for 
calculating rates to be compared with experiment. 

An example of the structure in the unconvolved cross 
sections is shown in Fig. [T] There we show raw results of 
the rovibration-only frame transformation calculation for 
J^=2, both using the full rovibrational Hamiltonian to 
calculate the rovibrational states, and using a rigid rotor 
approximation for the rovibrational states. The results 
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FIG. 7: Raw output of frame transformation calculation. 
Plotted is the dissociative recombination probability for the 
J"'"=2 vibrational ground state, s-wave channel. Top, result 
with full rovibrational Hamiltonian; bottom, result with rigid 
rotor states. 



are markedly different, showing that the strong mixing 
of rotation and vibration in LiH^, even at low J, affects 
the structure in the cross sections for individual entrance 
and exit channels. 

The first excited rovibrational state lies at 7.3meV. 
The sixth and ninth excited state, corresponding to ex- 
citation in the dissociative R direction and excitation in 
the 7 direction - rotation of the H2 - lie at 53meV and 
77meV, respectively. As is clear from Figure [71 there is a 
prominent series of narrow rydberg resonances converg- 
ing to the 53meV threshold, which serve to enhance the 
DR rate. It is therefore clear that excitation in the dis- 
sociative direction plays the largest role in the indirect 
DR process for this molecule, as opposed to rotational 
excitation or excitation in the H2 stretch coordinate. 

For the purpose of calculating rates to be compared 
with experiment, we find that the rigid rotor treatment 
is probably sufficient, though it apparently overestimates 
the cross section slightly. Not including the rotational 
frame transformation of Child and Jungen, we have com- 






cc 

Q 



10000 



1000 



100 



10 



Rotational + Rovibrational 
Rovibrational 
Rigid rotor 




0.0001 0.001 0.01 0.1 

Incident electron energy (eV) 



FIG. 8: Calculated dissociative recombination rate at 300°K, 
assuming experimental resolution of \/2 meV in the parallel 
and transverse directions. The rotational plus rovibrational 
curve is our final result, and the other two curves are from 
the rovibrational transform only, using rigid rotor or full rovi- 
brational states. 



pared the rovibration-only frame transformation using 
the full rovibrational states to that using the rigid rotor 
states. We find that the rigid-rotor treatment yields a 
DR rate consistently about 20% higher than the full rovi- 
brational treatment. The full calculation, employing the 
rotation/rovibrational frame transformation and the full 
rovibrational states, was not completed, due to numeri- 
cal difficulty. Instead, we perform the rotational trans- 
formation of Child and Jungen with the rovibrationally 
transformed S-matrix calculated from rigid-rotor states. 
On the basis of the comparison between the calculations 
not including the rotational transformation of Child and 
Jungen, we estimate that this treatment probably over- 
estimates the cross section by about 20%. 

Our convolved results are shown in Figure [51 Wc show 
the DR rate calculated at 300°K and including states up 
to J~^=9 (for the rovibration-only transformations) or 
J=ll (for the rotational and rovibrational transforma- 
tion). We show three results: from using the full rovibra- 
tional Hamiltonian, with no rotational transform; from 
using a rigid rotor approximation, with no rotational 
transform; and from using a rigid rotor approximation, 
with the rotational transformation of Child and Jungen. 
The former two calculations demonstrate the effect of in- 
cluding the full rovibrational dynamics, and as mentioned 
immediately above, the rigid rotor result exceeds the full 
rovibrational result by approximately 20 percent, which 
factor is fairly independent of the incident electron en- 
ergy. The latter calculation should be considered our fi- 
nal result, with the caveat that it probably overestimates 
the rate by about 20%. Nuclear statistics are included 
for the full rotational/rovibrational transformation, but 



not for the rovibration-only transformation, because the 
rovibration-only frame transformation destroys the per- 
muation symmetry of the overall wavefunction. The rates 
are comperable but a bit higher than the corresponding 
rates for H^ , by a factor of two or three. The effect of 
including the rotational part of the transformation is to 
further lower the results by about 10% in the low-energy 
region, and 50% in the high-energy region. 



VIII. CONCLUSION 

We have applied the method of Ref. [l| to the calcu- 
lation of the indirect DR rate of LiH^ + e~ . A cen- 
tral aim of our treatment was to analyze the effect of 
including the full rovibrational dynamics of the cation. 
We have found that although the full rovibrational treat- 
ment produces channel energies and unconvolved cross 
sections considerably different from a rigid rotor treat- 
ment, a rigid rotor treatment is amenable to the calcu- 
lation of convolved cross sections to be compared with 
experiment, although it probably overestimates the DR 
rate for a floppy molecule such as LiH^ by a small and 
energy-independent amount. 

The main approximation in the present treatment is 
the use of energy-independent quantum defects. For the 
calculation of indirect DR rates for the present system, 
this approximation is expected to be very good, because 
the amount of energy transferred from electronic to nu- 
clear motion is rather small due to the small dissociation 
energy of LiH J . Methods to accurately treat the energy 
dependence of the fixed-nuclei s-matrix within a frame 
transformation exist [671 . [68l . [69[ and may be applied to 
this system in future work. 

The calculations presented here demonstrate that the 
indirect mechanism provides a powerful mechanism for 



dissociative recombination of LiHo 



Future work 



will seek to analyze the branching ratios for two- and 
three-body dissociation and to further study the nature 
of the indirect DR mechanism. 
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